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The density-matrix renormalization-group method is used to study the electronic states of a two- 
chain Hubbard model for CuO double chains of PrBa2Cu40s. We show that the model at quarter 
filling has the charge ordered phases with stripe-type and in-line-type patterns in the parameter 
space, and in-between, there appears a wide region of vanishing charge gap; the latter phase is 
characteristic of either Tomonaga-Luttinger liquid or a metallic state with a spin gap. We argue 
that the low-energy electronic state of the CuO double chains of PrBa2Cu40s should be in the 
metallic state with a possibly small spin gap. 
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I. INTRODUCTION 



Geometrical frustration in the charge degrees of free- 
dom of the CuO double chains of PrBa2Cu40s has pro- 
vided an interesting opportunity for studying anomalous 
metallic states realized in quasi-one-dimensional (ID) 
electronic conductors. Experimentally, it has been re- 
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that the system shows metallic conductivity 



along the chain direction down to 2K and that it shows 
insulating conductivity along the directions perpendicu- 
lar to the chain above ~140 K. The system may thus 
be regarded as a quasi- ID correlated conductor except 
at low temperatures. Below ^140 K, the conductivity 
of the system shows a 2D behavior within the crystallo- 
graphic a-b plane, 2 indicating that the coupling between 
the chains is not negligible. 

Anomalous behaviors in PrBa2Cu40s have been ob- 
served in particular in its charge degrees of free- 
dom. A nuclcar-quadrupole-resonance (NQR) experi- 
ment^ showed that the temperature dependence of the 
spin-lattice (1/Ti) and spin-spin (I/T2) relaxations is 
anomalous, indicating the presence of very slow fluctua- 
tions of charge carriers. A polarized optical spectrum for 
E || a (electric field E parallel to the a axis) demon- 
strated that the Cu02 planes of this material main- 
tain the charge-transfer (CT) insulating state with a CT 
gap of 1.4 eV^ which is typical of that of cuprates, 
but the spectrum for E || b (chain direction) exhibits a 
broad peak structure at lo > 40 meV with a very small 
but nonzero Drude weight (~2% of the total weight). 5 
Tomonaga-Luttinger liquid (TLL) behavior was also sug- 
gested; from the power-law dependence dominating in 
the higher-energy part of the optical conductivity, the 
exponent K p was estimated to be ~ 0.24. 5 Moreover, 
an angle-resolved photoemission spectroscopy (ARPES) 
study of Zn-doped PrBa2Cu40s showed a power-law be- 
havior of the spectral function near the Fermi level, which 
is indicative of a TLL state, and a value K p ~ 0.24 was 



estimatedi^ Such a small value of K p suggests that the 
long-range Coulomb interaction between charge carriers 
should be very strong. It was also reported that the 
system of the CuO double chains is near quarter filling 
of holes. 5 ' 7 A good opportunity has thus been provided 
for studying the anomalous charge dynamics of strongly 
correlated ID electron systems such as charge ordering 
(CO), charge fluctuation, and charge frustration. 

A few theoretical studies have so far been made on 
the CuO double chains of PrBa2Cu40s. Using Lanc- 
zos exact-diagonalization technique on small clusters, Seo 
and Ogata£ studied similar models and proposed possi- 
bility of the appearance of a metallic phase due to ge- 
ometrical frustration; the system sizes used were how- 
ever very small, so that whether the first charge exci- 
tation is gapped or gapless in the thermodynamic limit 
was not clear. Zhuravlev et ali^ investigated a spin- 
less Fermion model at half filling, which corresponds to 
a quarter-filled extended Hubbard model in the limit of 
large on-site repulsion U, and found that a metallic state 
is realized due to the effect of frustration among Fcrmions 
when long-range repulsive interactions among them com- 
pete. Amasaki et aZ.ii calculated the optical conductiv- 
ity spectrum in the mean-field theory and showed that 
it agrees well with experiment only when the stripe-type 
CO is present; they thereby argued that the observed 
anomalous optical response may be due to the presence 
of the stripe-type fluctuations of charge carriers in the 
CuO double chains. 

Motivated by such development in the field, we here 
investigate the ground state of the extended Hubbard 
model at quarter filling, an effective model for the 
CuO double chain of PrBa2Cu40s- An isolated double- 
chain system is considered in this paper because to 
clarify its nature is essential before the 'dimensional 
crossover' is taken into account. We use the density- 
matrix renormalization-group (DMRG) method^ for cal- 
culating the ground and low-energy excited states of the 
model, as well as its charge and spin correlation func- 
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tions. 

We will show that our model in its ground state has 
the CO phases with stripe-type and in-line-type patterns 
in the parameter space, and in-between, there appears a 
wide region of vanishing charge gap; the latter phase is 
characteristic of either a TLL or a metallic state with a 
spin gap. The spin gap remains open unless the filling 
deviates largely from a quarter. We will then consider 
possible experimental relevance of our results. We will in 
particular argue that the low-energy electronic state of 
the CuO double chains of PrBa 2 Cu40g shou 
metallic state with a possibly small spin gap 

This paper is organized as follows. In Sec. ' 
duce a single-band Hubbard model defined on a ID lat- 
tice for the CuO double chain of PrBa2Cu40s and give 
some numerical details of the DMRG method used. In 
Sec. Ill, we present calculated results for the charge gap, 
phase diagram, spin correlations, charge correlations, and 
Luttinger parameters, and discuss anomalous properties 
of the liquid phase of our system. In Sec. IV, we sum- 
marize the results obtained and discuss their possible ex- 
perimental relevance. 




FIG. 1: (a) Schematic representation of the lattice struc- 
ture of the CuO double chain. The solid circles represent Cu 
(d x 2_ y 2) orbitals and the open circles represent O (p CT ) or- 
bitals. (b) The corresponding two-chain Hubbard model, (c) 
Geometrically equivalent single-chain Hubbard model. 



II. MODEL AND METHOD 

The structural element of PrBa2Cu40s that we con- 
sider is the CuO double chain with unit cell Cu204ji*£ 
the geometry of which is illustrated in Fig.^a). The rel- 
evant orbitals on the Cu and O sites are the 3d x 2_ y 2 and 
2p a orbitals, respectively. From the antibonding band of 
these two orbitals, we may extract an effective model, i.e., 
a single-band extended Hubbard model for low-energy 
electronic states of the CuO double chains (see Fig.EJb)), 
where the O ions are no longer explicitly present. Such re- 
duction of dimensionality allows us to carry out detailed 
numerical studies. One should also note that this two- 
chain Hubbard model is topologically equivalent to the 
single-chain Hubbard model with next-nearest-neighbor 
interactions as shown in Fig.^c). 

The model we consider is thus described by the Hamil- 
tonian 

H = - h ^(ct ffCj+ i CT + 4 +la c ia ) 

ia 

- t 2 ^{c\ a C l+2<J + 4 +2(J C la ) 



+ U ^n^mi 

i 




where we use the notations for the ID chain defined in 
Fig.njc). c\ a (c ia ) is the creation (annihilation) operator 

of a hole with spin a (=t, |) at site i and = <y itJ Ci a is 



the number operator. t\ and Vi are the nearest-neighbor 
hopping integral and Coulombic repulsion, respectively, 
which are for the zigzag bonds connecting two chains of 
the two-chain model (see Fig. n^b)). t% and V2 are those 
for the next-nearest-neighbor bonds, which are for the 
chain bonds of the two-chain model. We assume t\ -C £2 
because of the bonding angle minimizing the overlap the 
Cu id x 2_ y 2 and O 2p a wave functions for the hopping pa- 
rameter ti; in the following, we use the value ii/^2 = 0.1 
unless otherwise stated, by confirming that the sign of 
t\ does not change the results. We assume a value of 
the on-site Coulomb interaction U as U jti = 20, which 
is somewhat larger than a typical value for cuprate ma- 
terials; by using this value, however, we can investigate 
interesting features of our model with a wide range of V\ 
and V2. We believe that this choice does not change the 
essential features of our results. The filling of holes in 
the CuO double chains is reported to be n ~ 0.5^; we 
therefore restrict ourselves to the case at quarter filling 
unless otherwise indicated, i.e., n = \ corresponding to 
a density of one hole per two Cu sites. Note that the 
repulsions V\ and V2 are frustrating interactions for the 
charge degrees of freedom of the system in the sense that 
two holes can choose to sit on two of the three sites of 
a triangle of the lattice in mutually competing energy 
scales V\ and V2 (see Fig.^b)). If we assume the inter- 
site Coulomb repulsions to be inversely proportional to 
the intersite distance, we have the relation V\ = V2V2, 
which is not far from the line V\ = IV2 where the two 
repulsions V\ and V2 exactly compete at ti — t 2 = 0. 

We use the DMRG method?^ a powerful numerical 
technique for a variety of ID systems, whereby we can 
obtain very accurate ground-state energies and expecta- 
tion values for very large finite-size systems. We study 
the chains with up to 128 sites. The open-end boundary 
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conditions are used; when we calculate local quantities, 
they are actually the ones averaged over the chain and 
when we show correlation functions, they are the ones 
measured from the midpoint of the chain. We use up to 
rn = 2000 density matrix eigenstates to build the DMRG 
basis. One needs to keep relatively large values of m 
to obtain sufficient accuracies in our system because the 
long-range hopping term increases the truncation error 
rapidly; this is in particular the case when t\ ~ 0. Ex- 
trapolating the DMRG ground-state energies calculated 
for different m values to an energy for vanishing discarded 
weight^ we obtain the ground state energy and excita- 
tion gap which are accurate to parts in 10 _3 i2- Although 
the largest source of errors in our calculations are finite 
size effects, we are able to extrapolate a number of gap 
values measured for different system sizes to a value for 
an infinite system-size chain. 

In the following, we first present calculated results for 
the charge gap in order to show that the geometrical frus- 
tration indeed induces the metal-insulator (MI) phase 
transition; the ground-state phase diagram is shown in 
the parameter space of the inter-site Coulomb interac- 
tions. We then present calculated results for the den- 
sity distributions, spin-spin correlation functions, charge 
velocities, and the density-density correlation functions, 
in order to discuss possible anomalous behaviors of the 
metallic phase. 
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III. RESULTS OF CALCULATION 
A. Charge gap 

The charge gap of an L-site chain is defined by 

A C (L) = ^[E (N^ + l,N l + l,L) 

+ E (N^ -l,Ni- 1, L) - 2E (N h N u L)] (2) 

where Eo(Nf,Ni,L) denotes the ground-state energy of 
a chain of length L with iVj spin-up holes and JVj. spin- 
down holes. In Fig. |2Ja), we show calculated results for 
A c (L)/t 2 at Vi/t 2 = 5 for several V 2 /t 2 values and sys- 
tem sizes L up to 80 sites. We find that A C {L) decreases 
monotonically with increasing L, so that we can extrap- 
olate A C (L) to the thermodynamic limit systematically 
by performing the least-square fit of A C (L) to a second- 
order polynomial in 1/L. The extrapolated results A c /t 2 
(= A C (L — > oo)/t 2 ) are shown in Fig.^b). It is evident 
that for both small V 2 /t 2 1-56 and large V 2 /t2 ^ 5.95, 
A c /t 2 is finite, i.e., the system is insulating, and in a wide 
range of V 2 /t 2 , i.e., 1.56 < V 2 /t 2 < 5.95, A c /t 2 vanishes 
within the accuracy of the extrapolation. We find similar 
situations for other values of V\ /t 2 as well, which demon- 
strates that a stable metallic phase indeed exists between 
two insulating phases. 

When t\ — t 2 — 0, we readily find from a total-energy 
calculation that there are two ordered phases, of which 



FIG. 2: Charge gap in the extended Hubbard model [Eq. 
at quarter filling for Vi/t2 = 5 and t\jti = 0.1. (a) A c (L)/t2 
calculated as a function of the inverse system size 1/L for 
= (filled circles), 1 (open circles), 2 (filled triangles), 
6 (open triangles), 8 (filled squares), and 10 (open squares). 
Dashed lines are the quadratic fits in 1/L. (b) V2/t2 depen- 
dence of A c /t2 in the thermodynamic limit L — + 00. The cir- 
cles show the results obtained by extrapolation of the DMRG 
data in (a). The dashed lines are a guide to the eye. 



the phase boundary locates on a V\ — 2V 2 line. The 
phase in the V± > 2V 2 region is characterized by an in-line 
CO, i.e., all the holes locate on one of the two chains of 
the two-chain model, whereas the phase in the V\ < 2V 2 
region is characterized by a stripe-type CO, i.e., holes lo- 
cate alternately on sites of both of the two chains form- 
ing a stripe-like pattern along the chain direction. The 
stripe-type CO corresponds to the 2Uf charge-density- 
wave (CDW) state and the in-line CO corresponds to 
the 4/cp-CDW in the single-chain notation of Fig^c). 
Note that, if we assume the intersite Coulomb repulsion 
to be inversely proportional to the intersite distance, we 
have V\ = \[2V 2 and hence the stripe-type CO should be 
realized, but if we assume a faster decay of the intersite 
repulsion due to the effect of screening, the energies of 
the two phases can be more comparable. 

By making t\ and t 2 finite to introduce quantum fluc- 
tuations, there appears the metallic state between these 
two insulating phases, where the two intersite repulsions 




FIG. 3: Ground-state phase diagram of the model [Eq. Q] 
in the parameter space of the intersite Coulomb interactions. 
The values t\ — Q.U2 and U — 20i2 are assumed. The open 
(solid) circles indicate the phase boundary between the 2&f- 
CDW (4/cf-CDW) and metallic states, at which the calcu- 
lated charge gap A c vanishes. The solid line corresponds to 
V2/t2 = Vi/(2i2)- Details of the metallic phase are discussed 
in Sec. Ill E. 



V\ and V 2 are mutually competing, as we have shown 
above. The ground-state phase diagram of our model 
(see Sec. Ill E for detailed discussion) therefore includes 
three phases in the parameter space as shown in Fig. [3J 
(i) the insulating phase with the 4fcp-CDW (or in-line 
CO) realized when V\ is large, (ii) the insulating phase 
with the 2/c^-CDW (or stripe- type CO) realized when V2 
is large, and (iii) the metallic phase in-between. These 
three phases are discussed further in the following sub- 
sections. 
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B. Charge ordering 

Let us first consider the three phases from the view- 
point of ordering of charge degrees of freedom. We first of 
all note that the CO is a state with a broken translational 
symmetry; there are two degenerate ground states, and 
in the open-end boundary conditions we use here, one 
of the two ground states is picked out by initial condi- 
tions of the calculation. The CO state is thus directly 
observable in our calculations. 

For small values of V 2 /t 2 (< V 1 /2t 2 - 0(h)), the ef- 
fect of Vi/fa dominates and leads to a CO with wave 
vector q = Akp in the single-chain notation (which corre- 
sponds to the in-line-type CO in the two-chain notation) ; 
namely, the ground state of the system is dominated by 
configurations such as 

...01010101010..., 

where "1" stands for a site supporting one hole and "0" 
stands for an empty site. Fig. 0fa) shows an example of 
the Friedel oscillation appeared in the density distribu- 
tion (m) in the finite-size system of L = 128. One can see 
the 4fci?-CDW behavior clearly. The oscillations persist 
at the center of the chain and the difference between the 



FIG. 4: Calculated hole densities at site i for (a) the 4/cf- 
CDW (or in-line-type CO) at Va/<2 =0 (crosses) and 1 (cir- 
cles), (b) 2fc F -CDW (or stripe-type CO) at Va/ta = 6 (circles) 
and 8 (crosses), and (c) uniform metallic phase at Va/ta = 2. 
We assume Vi/t2 =5 and use the system of L = 128. Results 
for central 40 sites are shown in (a) and (b). 



results for V 2 /t 2 = and 1 is hardly visible. We calculate 
(rij) for several chain lengths and confirm that the system 
size with L = 128 is sufficiently large and the finite-size 
effect is negligible. We also calculate the density-density 
correlation function N(l) defined as 



N(l) = (n i+ ini) - (n l+l ) (n,-) 



(3) 



with n$ = mi + riii and find that N(l) tends to a nonzero 
constant value as the intersite distance I increases: e.g., 
N(l 00) ~ 10~ 7 at Vi/h = 5 and V 2 /t 2 = 1. This in- 
dicates the presence of the long-range ordered Akp-CDW 
state. 

At V 2 /t 2 = where the Akp-CDW is stabilized, most 
of the charges are located on one of the two chains of the 
two-chain model although there are some charges on the 
other chain due to quantum fluctuation by a small value 
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of ti/ti) at tt/tz = 0, the charges are disproportionated 
on one of the chains completely. With increasing V 2 /t 2 , 
the charge gap decreases rapidly (see Fig. EJb)), but the 
4fej?-CDW correlation is hardly suppressed until reach- 
ing a sharp decline occurred around the MI transition 
point. This transition is of the second-order due to the 
presence of a small value of t\\ if t\ = 0, this transition 
should be of the first-order. These results suggest that at 
V 2 /t 2 < 1.56 and V\jt% = 5 the system may be regarded 
as a half-filled extended Hubbard model since most of the 
charges are disproportionated. We note that, in the case 
of t\ = 1 and t 2 = 0, the suppression of 4fci?-CDW cor- 
relation by V 2 has been suggested in the weak-coupling 
renormalization-group study^ 

For large values of V 2 /t 2 (> Vi/2t 2 + 0(h)), V 2 /t 2 
enhances the 2fc^-CDW fluctuations and induces a CDW 
with wave vector q = 2k f (or a stripe- type CO in the two- 
chain notation). The ground-state wave functions have 
dominant configurations of the type 

...011001100110... 

where we use the same notations as above. Fig. 01b) 
shows an example of the Friedel oscillation appeared in 
the density distribution (m) in the finite-size system of 
L = 128. The 2fc^-CDW behavior is clearly seen. The 
calculated values of N (I) tends to a nonzero constant as 
the intersite distance I increases: e.g., N(l — > oo) ~ 10~ 8 
at Vi/t 2 = 5 and V 2 /t 2 = 8. This indicates the presence 
of long-range ordered 2fc^-CDW state. 

When V 2 /t 2 is much larger than Vi/t 2 , 2fc F -CDW (or 
stripe-type CO) is stabilized, i.e., in the single-chain no- 
tation, two occupied sites and two empty sites come al- 
ternately. As shown in the next subsection, two spins 
on the neighboring two sites form a spin-singlet pair. As 
V 2 /t 2 decreases, A c decreases in proportion to V 2 /t 2 , and 
its behavior is similar to that of the model without hy- 
bridization (ti = t 2 = 0); i.e., A c = 2(V 2 - V 2 2kF ) 
around the MI transition point V 2 c,2kF . This suggests 
that the charges disproportionate almost completely as 
soon as the charge gap opens. It can also be verified by 
slight reduction of the amplitude of the 2kp -CDW oscilla- 
tion that occurs when V 2 /t 2 decreases. This means that 
the effect of V\/t 2 directly weakens the effect of V 2 /t 2 
and melts the 2/cf-CDW state. This situation is similar 
to the MI transition of the quarter-filled extended Hub- 
bard model. The transition is of the second order for 
t 2 > 0; it is of the first order only at t 2 = 0. 

For intermediate values of V 2 /t 2 (or around V± ~ 2V 2 ), 
the system is in the metallic state. As expected, the am- 
plitude of the Friedel oscillation is quite small as shown in 
Fig. 0Jc) , which should vanish at the center of the chain 
when L — - > oo. In this phase, competition of two incom- 
patible fluctuations, i.e., the 2kp and AkF (or stripe-type 
and in-line-type) charge fluctuations, prohibit the system 
from forming CO states. However, it is useful to know 
what is the dominant correlation; here, we examine the 
long distance behavior of the density-density correlations 
N(l), in particular its period of oscillation. We thus find 




v 2 /t 2 

FIG. 5: Order parameters Xik F and X2k F as a function of 
V2/t2 for various system sizes at Vi/ta — 1- The filled (open) 
symbols are the results for xik F (X2fc F )- The up-triangles, 
circles, and down-triangles denote the results for L = 16, 48, 
and 80, respectively. 



that the 2fcp CO correlation is dominant for 2V 2 > Vi 
and the 4/cf CO correlation is dominant for 2V 2 ^ Vi , as 
will be discussed further in Sec. Ill D. 

In order to explain relationship between the MI tran- 
sition and CO formation, we introduce the order param- 
eters defined as 

XA kP = ^- E &lAkFTin i = t- E(- 1 ) in j w 

for the 4/c F -CDW and 

X2kF = ^-^ 2kFr ^ = z^E cos (^) ^ ( 5 ) 

for the 2A;f-CDW. Here, the sum is taken over central 
L e ff (= L — a) sites excluding a/2 sites at both ends 
of the system in order to eliminate edge effects; here we 
take a ~ 8 - 12. 

The obtained results at Vi /t 2 = 5 as a function of V 2 jt 2 
are shown in Fig.[S]for a number of system sizes. We find 
that as the system size increases the slopes of the curves 
around the critical points become steeper and that the 
result at L — 48 is nearly equal to the result at L = 80. 
Thus, the situation in the thermodynamic limit may well 
be expected except near the critical point where the size 
effect is not negligible. Let us focus on the 2fcp-CDW 
transition at V 2 /t 2 ~ 5.95. Then, X2k F shows a sharp 
development around V 2 /t 2 w 5, which is rather smaller 
than the value of V 2 ' 2kF /t 2 shown in Fig.[5Jb). However, 
in thermodynamic limit, this turning point should agree 
with V 2 ' 2kF /t 2 , and the transition seems to be still con- 
tinuous. We note that the transition becomes sharper for 
larger Vi/t 2 (or smaller £2) values and is of the first-order 
in the limit V 2 — > 00 (or t 2 — > 0). A similar interpretation 
can also be given to the 4fcj?-CDW transition. 
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FIG. 6: (a) Spin-spin correlation function S(l) as defined in 
Eq. © at V 2 /t 2 = (crosses) and 1 (circles) with Vi/tz = 5. 
The distance I is measured along the ti-chain in units of the 
interval between two sites. The log-log plot in the inset shows 
that the power-law decay of correlations, (b) S(l) at Vilti = 8 
(circles) and 10 (crosses) with Vxjti = 5. The semilog plot in 
the inset shows the exponential decay of correlations. 



C. Spin correlations 

In order to study the spin degrees of freedom of the sys- 
tem, we now calculate the spin-spin correlation function 
defined by 

S(l) = (SfSf +l ), (6) 

with Sf = (rij-f — riijJ/2. 

Let us first discuss the 4fc_p-CDW phase (or in-line CO 
phase in the two-chain notation). Fig. |Sfa) shows the 
spin-spin correlation functions calculated at V 2 /t 2 = 
and 1 with Vi/t2 = 5 for the system of L — 128. The site 
i in Eq. |JB} is fixed at the middle of the chain. One can 
clearly find the antiferromagnetic correlations between 
neighboring spins localized on one of the two-chains of 
the two-chain model, i.e., formation of the 2kp spin- 
density- wave (SDW). One also notes that, because the 
charge fluctuations through the zigzag bonds connecting 
two chains are nearly zero due to a very small value of 
ti/t2, no spin correlations between two chains of the two- 
chain model are seen. As shown in the inset of Fig. EUa) , 
the spin-spin correlations decay as a power law, similar 
to those for the ID Heisenberg chain, which depends very 
little on V2 /*2 unless the 4fcp-CDW state is broken. This 




FIG. 7: Spin gap A s /t2 as a function of t\/t2 for V2/t2 = 
(crosses), 5 (squares), and 8 (open circles) at Vi/t<2 = 5 in 
the thermodynamic limit L — » 00. The dashed line is A s = 

0. 5(^A 2 ) 2 

suggests that the system may be well described as the 
ID Heisenberg chain combined with an empty chain, as 
the in-line CO state may suggest. 

Now let us turn to the correlations in the 2/c^-CDW 
phase (or the stripe- type CO phase in the two-chain nota- 
tion). Fig.|ijb) show the spin-spin correlation functions 
for V 2 /t 2 = 8 and 10 at V x jt 2 = 5. The behavior of S(l) 
for / = An (n — 0, 1, 2, ■ • ■) shows the antiferromagnetic 
correlations. This corresponds to the 2fcp-SDW state of 
one of the two chains at quarter filling. We find simulta- 
neously that there appears an SDW on the other chain; 
the phase of the SDW is shifted by 3tt/A from the other 
SDW. It might be as if the two SDW states coexist. 

Because a sign of S(4n) is opposite to those of S(An ± 
1), the spins on the neighboring two sites should be oppo- 
site. The semilog plot in the inset of Fig. OJb) shows the 
exponential decay of the spin-spin correlations. The cor- 
relation lengths are £ ~ 4.17 and 3.57 for ^2/^2 = 8 and 
10, respectively. S(4n — 1) and S(4n — 3), which corre- 
spond to the spin-spin correlations between two i2-chains, 
decay with the same correlation length as 5(4n). These 
behavior are similar to the case of the two-leg Heisenberg 
(and half-filled Hubbard) ladder system. The correlation 
length is close to that for the two-leg isotropic Heisenberg 
system £ - 3.19i^ This means the effective hopping be- 
tween the next-nearest-neighbor sites along a i2-chain, 

1. e., sites i and i + 4, is of the order of t\. Furthermore, 
we have found that S(l) decays exponentially even in 
metallic regime near the 4fcp-CDW insulating phase. 

From the exponential decay of the spin-spin correla- 
tions, we expect the system to be in a gapped spin-liquid 
state. In order to confirm it, we have calculated the spin 
gap defined by 

A s (L) = E (Ni + l 1 N l -l,L)-E a {N^-l,N l + l,L) (7) 

for an L-site chain. In Fig. [7| we show the extrapolated 
values A s (= A S (L -> 00)) at Vi/t 2 = 5 for V 2 /t 2 = 
(4£: F -CDW phase), 5 (metallic phase), and 8 (2/c F -CDW 
phase). In the 4&F-CDW phase, we find that A s is al- 
ways zero, which is quite natural since the system can 
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FIG. 8: Density-density correlation function N(l) defined in 
Eq. plotted on a log- log scale. The results are for V>2/t2 ~ 3 
(filled triangles), 4 (open squares), and 6 (filled circles) at 
Vi/t2 = 5, with L = 128. The distance I is measured along 
the ID chain of Fig. ^c) in units of the interval between two 
sites. 



TABLE I: Calculated values of the Luttinger parameter K p 
and charge velocity v c at Vi jti =5. 



V a /t2 


K P 




2 


0.55 


0.295 


3 


1.25 


0.738 


4 


0.88 


0.401 


5 


0.27 


< 0.01 



even for finite interactions, one may expect that the sys- 
tem at low energies should behave as a two-band models 
In the two-band model, it is known that only 4kp-CDW 
and d-wave superconducting (SCd) correlations decay as 
power laws, and all other correlations decay exponen- 
tially; i.e., either SCd correlation or 4fci?-CDW correla- 
tion dominates for both weak and strong interactions. 18 
The 4/cp-CDW correlation and SCd correlation at large 
distances can be written like 

N(l) «cos(4M)r 2K <>, (11) 



be presumed to be a ID Heisenberg model as mentioned 
above. In the 2fcj?-CDW phase, we find that the spin gap 
is always finite except at t\ — and the size of the gap 
can be scaled with t\, i.e., 



A, 
ta 



2 U 2 



(8) 



as expected from the assumption that the neighboring 
two sites couple with the exchange interaction J. Because 
the spin gap for a two-leg isotropic Heisenberg ladder 
is A s ~ J/2j— we can estimate the effective exchange 
interaction as 



J 



(9) 



As for the metallic phase, A s goes through a maximum at 
a finite value of because an increase of ii/i2 leads to 
an increase in the effective exchange coupling, as well as 
a decrease in the 2/cf-CDW fluctuations. This behavior 
reminds us of the spin gap of the two-leg CuO ladder 
materials where the gap behaves similarly as a function 
of the energy difference between the O and Cu sites^ 



D. Luttinger parameter 

Let us calculate the Luttinger parameter K p which en- 
ables one to know competitions among correlations in the 
metallic phase. The Hamiltonian Eq. has the energy 
dispersion relation 



e(k) = —2ti cos k — 2t2 cos 2k. 
in the absence of interaction. If t\ <C £2 and 



(10) 
the 



12 anu n = \. 

Fermi surface has four Fermi points i/c^ and ±fcj? 2 , and 



and 



with 



(AL,A 



-1/2K„ 



Aj — -j= {^Ci+ii - Cjxc 1+ i T ) 



(12) 



(13) 



respectively. 

In Fig. |S1 we show the density-density correlation func- 
tion N(l) defined in Eq. Q on a log-log scale in order to 
estimate the exponent of the decay. At large distances 
(I > 10), the envelope of the correlation is consistent 
with a straight line on the log-log scale and therefore one 
may say that it decays as a power law. Fitting the de- 
cay to Eq. Ijlll) allows us to determine the value of K p 
although the fitting includes some errors due to Friedel 
oscillations in the data. The values obtained are given in 
Table which would be of sufficient accuracy when we 
consider available experimental results. 

At V2/t2 — 2.5, N(l) decays most rapidly and K p is the 
largest. This means that the long-range interactions are 
effectively weakened since V\ and V2 are canceled with 
each other. Such a situation seems to be realized around 
^2/^2 = Vi/(2t 2 )- Generally, long-range Coulomb re- 
pulsions suppress the value of K p . In our results, K p 
decreases when the values of V\ and V2 deviate from 
the relation V\ = 2V2, whereby the effective interaction 
strength increases. Now we note that the 2fcp-CDW tran- 
sition occurs when K p = 1/4, which may be understood 
if we suppose that two independent ^-chains at quarter 
filling undergo a '4fcf'-CDW transition. On the other 
hand, the 4fcp-CDW transition occurs when K p = 1/2. 
This may be understood if we assume that the system 
consists of a half-filled chain and an empty chain in the 
4fe F -CDW phase. 
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We note that the SCd correlations decay more rapidly 
than the 4fcp-CDW correlations in the whole metallic 
phase even if K p > 1. The inverse relationship of ex- 
ponents predicted by a bosonization pictur e) 19 ! 20 seems 
not to be preserved. The same situation was also seen 
for the two-leg Hubbard ladder^ where it is known that 
enhancement of the SCd correlation is not adequately 
observed due to the finite-size effect, i.e., the widely sep- 
arated energy levels near the Fermi energy. 

The slopes at 1/L -> in Fig. Ufa) yield nv c /(2K p ), 
from which we can estimate the charge velocity v c . We 
find that v c goes through a maximum around V 2 /t 2 = 
2.5 where the electrons can move almost freely despite 
the effect of double occupancy on a site. We note that 
the calculated very small charge velocities near the 2/cf- 
and 4&F-CDW states (see Table may suggest heavy- 
Fermionic features, for which the geometrical frustration 
may play a crucial rolem^MiSl 



E. Phase diagram 

In Fig. |3 we show the ground-state phase diagram 
of the system in the parameter space of the intersite 
Coulomb interactions Vi/t2 and V^/ta. The phases are 
distinguished by the presence of a gap for spin and/or 
charge excitations. Naturally, a CDW state is expected 
if cither V\ / t 2 or V 2 /t 2 is large. The ground state of the 
2fc F -CDW (4fc F -CDW) phase is a CO state with (with- 
out) spin gap and always insulating. Then, if both V~i/t 2 
and V 2 /f 2 are sufficiently small or the effective long-range 
Coulomb interaction is weak due to frustration between 
two interactions, the system may be in a metallic phase. 
The metallic phases can be divided into two regions. One 
is the spin-gapped-liquid (SGL) phase where there is a 
finite spin gap but no charge gap. Here, the system is 
under the influence of 2fcp-CDW transition and the 2kp- 
CDW fluctuation is dominant. Another is the TLL phase 
where neither spin gap nor charge gap is present and the 
4/c^-CDW fluctuation is dominant. 

For a lucid discussion, let us now assume the case 
U/t 2 — > 00 and consider some limiting cases. When 
V\/t 2 = 0, the 2fcp-CDW transition occurs at V%lt 2 — 
2.3. If the system can be regarded as two separate t 2 - 
chains due to very small ti/t 2 , this transition would be 
equivalent to the 4/cf-CDW transition of the quarter- 
filled spinless Hubbard chain with the nearest-neighbor 
hopping t 2 and the next-nearest-neighbor Coulomb re- 
pulsion V 2 , where the critical value of V 2 /t 2 is expected 
to be 2, which is consistent with our results. 

When V 2 /t 2 = 0, the Akp-CDW transition occurs 
around V 2 /t 2 ~ 2. We can estimate this value from 
comparison between the energies of the metallic and in- 
sulating solutions. First, the metallic solution may be 
obtained as follows. Consider the system consisting of 
two independent quarter-filled i2-chains. Then, if we as- 
sume the Coulomb repulsion between the i2-chains, the 



total ground-state energy may roughly be evaluated as 

2E%{n=\) + ^-L. (14) 

Here E 1 ^ (n = g) is the ground-state energy of a spinlcss- 
Fermion model defined on a single chain at quarter filling 
with the nearest-neighbor hopping t%, and Vi/4 is the av- 
erage Coulomb interaction. Next, the insulating solution 
may readily be obtained because it corresponds to the 
4fcp-CDW state, i.e., one of the ^-chains is at half fill- 
ing and another ^2-chain is empty, whereby we obtain 
the ground-state energy to be almost zero. Hence, the 
condition for emergence of the Akp-CDW state is that 
the value of Eq. I|14|l is larger than zero. If we use an 
estimation Eq 2 (7i = i) ~ —0.25t 2 L made from the Lieb- 
Wu equations^ the critical interaction strength can be 
obtained as V\jt 2 — 2. We should note that, because we 
have used a small but finite value of t\/t 2l the actual crit- 
ical values of the interactions calculated are rather larger 
than the above estimation. 

Next, we turn to the case where both Vi/t2 and V2/t 2 
are nonzero. Naturally, the system is metallic for small V± 
and V 2 (;$ 0(t 2 )), and the metallic phase becomes nar- 
rower as the long-range Coulomb interactions increase 
along the line V 2 /t 2 = Vi/(2t 2 ). For very large values 
of V\ and V 2 , the boundary between the 2fcp-CDW (the 
4&F-CDW) and the SGL (TLL) phases becomes parallel 
to the line ^2/^2 = Vi/(2t2), and the distance between 
the two phases remains to be of the order of t 2 (t±) even 
in the limit of Vi —* 00 and V 2 — * 00. Even if there is no 
hybridization, i.e., t\ = t 2 = 0, the charge mode is gap- 
less on the line V 2 /t 2 = V\/(2t 2 ). This means that the 
effects of V\ and V 2 are precisely cancelled out on that 
line. Now, if finite values of t\ and/or t 2 are introduced, 
the energy gain due to moving holes expands the width of 
the line. We point out that the width of metallic regime 
depends on the difference SV between two repulsive inter- 
actions V\ and V 2 (SV = V\ — V 2 ), rather than the values 
V\ and V 2 themselves. To check this notion, we calculate 
the critical strength of V 2 / 1 2 for very large value of Vi/t 2 . 
The critical interaction strengths for the 4/cf-CDW and 
2fci?-CDW transitions are obtained respectively as fol- 
lows: V 2 /t 2 fa 249.5 and 253.5 for V x /t 2 = 500, and 
^2 A2 ~ 499.5 and 503.5, where we assume V x /t 2 = 1000 
and U/t 2 = 4000. If h = (t 2 = 0), the boundary of the 
2k F -CBW (the 4fc F -CDW) phase should be asymptotic 
to the line V 2 /t 2 = V l /{2t 2 ). 



IV. SUMMARY AND DISCUSSION 

Using the DMRG method, we have studied the elec- 
tronic states of a two-chain Hubbard model at quarter 
filling that models the low-energy electronic states of the 
CuO double chains of PrBaaCu-jOs- We have shown that 
the model has the CO phases with the stripe-type and 
in-line-type patterns in the parameter space and between 
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FIG. 9: Normalized optical conductivity on a log-log scale. 
Parameter values used for calculation are: U/t<2 = 20, Vi jti = 
2, Vi/t2 = 5, L = 32, and n = 0.46875. Finite broadening 
width r\ = 0.4^2 is introduced. The solid line shows a fitting 
to the form <j(lu) oc aj~ 7 . 7 ~ 1.5 is found. 

the two CO phases there appears a wide region of vanish- 
ing charge gap, which is due to geometrical frustration 
of the long-range Coulomb interactions. The phase dia- 
gram of our model consists of four phases characteristic 
of (i) the stripe- type CO (or 2/cf-CDW), (ii) in-line-type 
CO (or 4fcF-CDW), (iii) Tomonaga-Luttinger liquid, and 
(iv) spin-gapped liquid. For the liquid phases, we have 
suggested possible emergence of a heavy- Fermion-like be- 
havior, which can occur near the region of the metal- 
insulator phase boundary. 

Finally, let us discuss consistency between our nu- 
merical results and available experimental data on 
PrBa2Cri408. Comparable experimental features are the 
appearance of the TLL behaviors and estimation of the 
value K p ~ 0.24i£i& We can explain such a small value 
of K p only by assuming that the system is in the 2fcp- 
CDW state or in the SGL state near the 2fcp-CDW phase 
boundary. In these two phases of our phase diagram, 
however, a small but finite spin gap clearly exists except 
at ii = 0. We have also confirmed that the spin gap 
remains open even if the system is doped slightly away 
from quarter filling. Thus, the phase is not TLL unless 
t\ = 0, which is apparently not consistent with exper- 
iment. One possibility for reconciliation may be that 
the CuO double chain system has a very small value of 



ti, so that experiments performed with finite resolution 
cannot detect the effects of the spin gap. Another pos- 
sibility is that the optical and ARPES spectra observed 
in relatively high energies are insensitive to the opening 
of the spin gap. To demonstrate this, we have calculated 
the optical conductivity of our model with the dynam- 
ical DMRG method^ and confirmed that the shape of 
the higher energy part {hu> > A c ) of the conductivity 
does not depend on the existence of small values of t\. 
In Fig. El we show an example of the ii-dependence of 
the optical conductivity calculated for a possible set of 
parameters. We find that the spectra for t\/t2 = and 
0.1 are almost the same and the power-law decay can be 
seen at Tito > A c . From a fit of the results to the form 
<j(uj) oc w -7 , we estimate the value 7 ~ 1.5, which cor- 
responds to K p ~ 0.22. This value is very close to the 
value obtained experimentally. Here we point out that 
the long-range CDW ground state cannot reproduce the 
experimental spectrum but a small hole doping of only a 
few % changes the shape of spectrum largely leading to 
the power-law decay of the spectra at high energies. We 
therefore suggest that the actual material is located in 
the SGL phase near the 2/cf-CDW phase or in the 2kp- 
CDW phase with slight doping. We also want to suggest 
that, if the NMR Knight-shift measurement is made on 
the spin degrees of freedom of the double-chain Cu(l) site 
of PrBa2Cu408, one should be able to detect anomalous 
behaviors associated with the spin-singlet formation due 
to the stripe-type fluctuation of charge carriers. 
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